How Sensitive are Di-Leptons from p Mesons to the High Baryon Density Region? 
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We show that the measurement of di-leptons might provide only a restricted view into the most 
dense stages of heavy ion reactions. Thus, possible studies of meson and baryon properties at 
high baryon densities, as e.g. done at GSI-HADES and envisioned for FAIR-CBM, might observe 
weaker effects than currently expected in certain approaches. We argue that the strong absorption 
of resonances in the high baryon density region of the heavy ion collision masks information from 
the early hot and dense phase due to a strong increase of the total decay width because of collisional 
broadening. To obtain additional information, we also compare the currently used approaches to 
extract di-leptons from transport simulations - i.e. shining, only vector mesons from final baryon 
resonance decays and instant emission of di-leptons and find a strong sensitivity on the method 
employed in particular at FAIR and SPS energies. It is shown explicitly that a restriction to p 
meson (and therefore di-lepton) production only in final state baryon resonance decays provide a 
strong bias towards rather low baryon densities. The results presented are obtained from UrQMD 
v2.3 calculations using the standard set-up. 
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Quantum-Chromo-Dynamics (QCD) predicts that the 
properties of hadrons change when they are brought into 
a (hot and/or dense) nuclear environment [l], @|. This 
modification is due to the interaction with the surround- 
ing medium which eventually leads to chiral symmetry 
restoration at high baryon densities and/or high temper- 
atures 0] . The experimental verification of this theoreti- 
cal prediction is one of the most challenging questions in 
modern strongly interacting matter physics. 

Among the non-strange mesons the p meson plays a 
dominant role in these investigations. It has a short life- 
time and therefore it has a large probability to decay 
inside the reaction zone when created in heavy ion colli- 
sions. It couples strongly to nuclear resonances and, most 
important, it has a non- negligible chance to decay into di- 
leptons which leave the interaction zone essentially with- 
out any further interaction. Thus, the di-lepton channel 
seems to offer a unique chance to study the high baryon 
density properties of the p meson. 

Theoretically the question of how the spectral func- 
tion of the p meson changes in the medium is still under 
active discussion. There is certain theoretical evidence 
that the p meson is broadened if put into the nuclear 
medium In contrast, Hatsuda and Lee predicted 

a lowering of the p meson mass in a nuclear environment 
based on QCD sum rules calculations A result which 
has also been found by Brown and Rho @, 0, HI- On 



the other hand, more recent calculations indicate that 
the pole mass of the p meson remains almost unchanged 
in the nuclear medium [l(| [H|. However, these cal- 
culations rely on specific assumptions on the coupling 
strength of the p meson to the nuclear resonances and on 
the branching ratios whose validity can presently only 
be proven by comparison to experimental data. For the 
present status of the theoretical spectral function calcu- 
lations for vector mesons we refer to [HI, [l]| [l4|, [l5|, . 

Experiments have been launched to verify these theo- 
retical predictions. In proton- nucleus collisions [13, EH 
at 12 GeV a decrease of the p meson mass with increasing 
baryon density ps as m(ps)/m(0) = 1 — 0.09ps/ pbo - 
about half of the value predicted by theory - but without 
an increase of the p meson's width has been reported. 
The CLAS collaboration reports that the experimental 
data of photon-induced reactions is compatible with no 
shift of the p meson pole mass and no additional broaden- 
ing to the theoretically estimated collisional broadening 
In contrast, the di-lepton data in In+In collisions 
at 158 AGeV [2(| are best described using essentially the 
free p meson pole mass but a considerable broadening of 
the spectral function. At lower energies of 2 AGeV (see 
[2H for the newest data at 1 GeV), the HADES collabora- 
tion has recently published di-lepton spectra [22|. Here, 
a deviation from the yield calculated from a hadronic 
cocktail fit in the vicinity of the p meson mass is visible 
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but due to the many sources of di-lcptons there is no con- 
clusive explanation for that suppression yet. How much 
of these different experimental findings can be exclusively 
attributed to the different environments, i.e. cold nuclear 
matter in proton-nucleus reactions, an expanding meson 
dominated fireball after a possible phase transition from 
a quark gluon plasma in high energy nucleus- nucleus col- 
lisions at the SPS or a baryon dominated expansion in 
reactions at about 2 AGcV is still a matter of debate. 

To link the final state di-lepton data to the in-medium 
spectral functions of the hadrons detailed quantitative 
theoretical simulations of the baryon density distribution 
at the p meson production and decay /absorption point 
are necessary This allows then to calculate di-lepton 
spectra from the simulations. Unfortunately, up to now 
different approaches are used to convert the calculated 
hadron spectra into di-leptons. In general four different 
approaches can be identified: 

• Explicit propagation of stable particles, baryon 
and meson resonances, decays of resonances into 
other mesons (especially baryonic resonances into p 
mesons), as well as irn — > p scattering. Di-lcptons 
are emitted continuously from vector mesons and 
baryon resonances with their respective locally 
given total width while the resonances are prop- 
agated ("shining") [11, [II 11, [II III, 

• Explicit propagation of stable particles, baryon 
and meson resonances, decay of resonances into 
other mesons (especially baryonic resonances into p 
mesons), as well as mr — > p scattering. Di-lcptons 
are only emitted at the point of decay (not absorp- 
tion) of the hadronic resonance [28[ , 

• Explicit propagation of nucleons, pions and nucleon 
resonances while vector meson degrees of freedom 
are not propagated explicitly Di-leptons are pro- 
duced via the eVMD model from nucleonic reso- 
nances at the point of the decay of the nucleonic 
resonances. This model assumes di-lepton produc- 
tion via intermediate p/ui states [Itl |29|. 

• Explicit propagation of nucleons, pions, kaons and 
A resonance. Di-lcptons are produced by decay 
of nucleonic resonances according to the branching 
ratios. This is supplemented by the di-lepton pro- 
duction from mesonic resonances, which are calcu- 
lated by folding the yfs distribution of the elemen- 
tary nucleon-nucleon collisions with the mesonic 
production cross sections and the corresponding 
branching ratios into di-leptons [3tjj . 

While the first method is sensitive to all stages of the 
collision and allows for a dynamical treatment of the col- 
lisional broadening, the other methods rely on the actual 
decay of the meson or baryon resonance. Thus, the dif- 
ferent approaches probe different baryon density regimes 
and might therefore provide different results for the ex- 
tracted di-lepton rates. Another problem is posed by 
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FIG. 1: (Color online) Time evolution of the local rest frame 
baryon density ps averaged over the positions of the individ- 
ual hadrons for central Au+Au/Pb+Pb reactions at various 
beam energies. 



the implementation of brcmsstrahlung especially at low 
beam energies. This discussion has recently been revived 
by the calculations of [H M EL H . 

In this situation it is necessary to study the general 
differences between the above discussed approaches and 
explore the baryon density probed by the p meson in the 
FAIR energy regime. This helps to provide a theoret- 
ical error margin for further detailed model studies on 
the change of the in-medium spectral functions at these 
energies. 

We perform this study for massive nuclear reactions in 
the energy range of 2A GeV < E lah < 30A GeV. This 
range marks the expected transition towards chiral sym- 
metry restoration, but also the transition from baryon 
dominated to meson dominated matter. Dedicated facili- 
ties to explore this energy domain are the FAIR project at 
GSI [H[ and the critRHIC program at BNL [pi. For our 
studies we employ the UrQMD(v2.3) model [IE 11, E}. 
It is a non-equilibrium transport approach based on the 
covariant propagation of hadrons and strings. All ele- 
mentary cross sections are calculated by the principle of 
detailed balance, from the additive quark model or are 
derived from available data. This model allows to study 
the full space time evolution for all hadrons, resonances 
and their decay products. It permits to explore the emis- 
sion patterns of the resonances in detail and to gai n in- 
sight into the origin of the resonances [1^, H, EL I42T [H[ . 

In this approach the resonances are treated with their 
vacuum properties during the whole evolution of the re- 
action. However, one should note that the particle prop- 
erties are dynamically modified in a hot and/or dense 
medium due to the coupling of the p meson to the sur- 
rounding hadrons (at SIS energies especially the baryon 
resonances are important [H[). The p meson is assumed 
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FIG. 2: (Color online) Gain and loss rates of p mesons 
separated gain from collisions ("gain(collisions)"), gain 
from decay ("gain(decays)") and loss due to absorption 
("loss(absorbed)") and loss in decays ("loss(decayed)"). From 
top to bottom we display central Au+Au/Pb+Pb collisions 
at 2, 11 and 30 AGeV. 



FIG. 3: (Color online) Gain and loss of p mesons separated for 
collisions and decay. From top to bottom we display central 
Au+Au/Pb+Pb collisions at 2, 11 and 30 AGeV. We show as 
well the difference of gain and loss, the number of p mesons 
present in the system as a function of time. From top to 
bottom we display central Au+Au/Pb+Pb collisions at 2, 11 
and 30 AGeV. 



to have a lifetime according to an exponential distribu- 
tion with a mean lifetime r of l/r po ; e ~ 1/150 MeV, in 
addition collisional broadening is implicitly taken into ac- 
count for the calculation of di-leptons by the interaction 
of the p meson with the evolving medium. 

UrQMD has been successfully applied to study light 
and heavy ion reactions from -Ebcam = 90^4 MeV to 
y/sNN = 200 GeV. The present version of the trans- 
port model is able to describe the yields and the trans- 
verse momentum spectra of different particles in proton- 
proton, proton-nucleus and nucleus-nucleus collisions 
[371 144| | . Detailed comparisons of UrQMD with a large 
body of experimental data at SIS and FAIR energies can 



be found in [37], IH[ • For further details of the model the 
reader is referred to [35|, HH, [37j • 

Let us start the discussion by displaying the time evo- 
lution of the baryon density in central Au+Au/Pb+Pb 
reactions at 2, 11 and 30A GeV (see Fig. [T]). The baryon 
density is averaged over all hadron positions and is cal- 
culated locally in the rest frame of the baryon current 
(Eckart frame) averaged over the position of every baryon 

as pb = j with j ^ = (pb , 0) ■ Details on the calcula- 
tion of the baryon density are discussed in the Appendix. 
Note that the maximal baryon density grows with in- 
creasing beam energy. The question to be asked is, how 
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FIG. 4: (Color online) Baryon density distribution at the 
space points where the p mesons decay. From top to bottom 
we display central Au+Au/Pb+Pb collisions at 2, 11 and 30 
AGeV 



sensitive are di-lepton obscrvables on the high baryon 
density stage of the collisions? 

In order to investigate this question, we first review 
the different production and loss mechanisms for the p 
meson. In UrQMD the p meson can be produced from 
the decay of a high mass (meson or baryon) resonance or 
directly in a collision of two particles (e.g. 7r + tt — > p) 
which includes also the production from string fragmen- 
tation. The p meson is destroyed by two different mech- 
anisms. It can decay ("loss (decayed)", e.g. p — > h + -k~) 
or it can be absorbed in collisions ("loss (absorbed)", e.g. 
p + 7r — > ai). 

In Fig. [2] we show the respective gain and loss rates 
of p mesons separated for the different contributions for 
central Au+Au and Pb+Pb reactions at = 2, 11 



FIG. 5: (Color online) Distribution of the baryon density at 
which the e + e~-pairs from the p vector meson are emitted. 
The thin lines denote calculation where the full branching 
ratio into di-leptons is attached to both the p decay and 
the p absorption vertices. The thick lines indicate the re- 
sults obtained by the shining method with continuous emis- 
sion of di-leptons. From top to bottom we display central 
Au+Au/Pb+Pb collisions at 2, 11 and 30 AGeV. 



and 30A GeV. The different processes [46[ in which a 
p meson is produced are denoted as "gain (collisions)" 
or "gain (decays)" (meaning the stem from the decay of 
another resonance) . The loss term differentiates also be- 
tween p mesons which have decayed "loss (decayed)" and 
those which are absorbed in collisions "loss (absorbed)" 
as discussed above. 

At 2A GeV, one observes that the production rate of p 
mesons is dominated by the decay of resonances (~ 80%, 
full line) as compared to the formation in a s-channel 7T7t 
scattering (~ 20 %, dotted line). This observation is in 
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line with the expectations that a dominant production 
channel for the p meson in low and intermediate energy 
heavy ion collisions is the decay of baryon resonances 
[38| . A previous detailed analysis of the p meson produc- 
tion channels at 2A GeV within the UrQMD approach, 
found in [43j | confirms this interpretation in detail. The 
maximum of the p meson production rate coincides with 
the maximum baryon density around t ~ 9 fm. How- 
ever, these p mesons are subject to frequent interactions 
with the surrounding baryons resulting in rather short 
life times of the p mesons as indicated by the large ab- 
sorption rate (dashed line). Only towards the end of 
the high density stage (when the p absorption processes, 
e.g. p + B — > B*' cease) p mesons can decay directly 
as denoted by the dashed-dotted line. In the present 
model p meson absorption accounts for the main loss of 
p mesons, while the decay accounts for only 30% of the p 
meson loss. It is clear that these features might lead to 
different time-dependent di-lepton yields, depending on 
the method with which the di-lepton rates are extracted 
from the numerical simulation. 

At higher energies (11 A GeV, 30,4 GeV) this low en- 
ergy line of arguments changes. Here one observes two 
distinct phases for the production and decay/ absorption 
of the p meson. Initially p meson production from colli- 
sions proceeding either via string formation and fragmen- 
tation or via meson-meson scattering dominates the gain 
term (dotted line). However, also the absorption proba- 
bility is rather high in this stage of the reaction resulting 
in a quick re-absorption of the p meson (dashed line) . For 
di-lepton calculations it becomes clear that only a "shin- 
ing" approach has the potential to provide information 
on this stage, whereas the approaches which depends on 
the decay of the resonance do not allow to extract this 
information. However, also at higher energies, the pro- 
duction of p mesons from resonance decays in the late 
stage of the reaction is sizeable. This is evident because 
the sequential processes NN —>B*+X,B*—>p + X 
will need a certain time and therefore triggers on later 
stages of the reaction. Thus, the p meson production 
from baryon resonance decays again leads towards a self- 
triggering of p meson production and subsequent decay 
at rather moderate densities, even at energies of 11 A GeV 
and 30A GeV. 

To gauge our assumptions for the absolute importance 
of the different processes the integrated p meson rates 
are shown in Fig. [3j The nomenclature is the same that 
has been used in Fig. [51 The temporal evolution of the 
integrals for three values of beam energies from £7i a b = 
2A GeV to £i a b = 30A GeV is calculated. In addition, 
the sum of the gain (S ga in) and the sum of the loss terms 
(Eioss) are displayed. If one subtracts the sum of the loss 
terms from the sum of the gain terms one gets the yield 
(n(t)) of p mesons which is present at each time step of 
the collision. 

As discussed above, at the lowest energy the gain via 
resonance decays dominates over the gain from direct p 
meson production due to kinematical constraints. How- 



4.5 

4.0 

3.5 


ftJec 

-- ^created 


3.0 




§2.5 

CD 

Q.2.0 


— 


1.5 




1.0 




0.5 


f 2 AGeV (Au+Au) 


4.0 

3.5 




H 1 — 1 — • — 1 — • — 1 — • — 1 — • — 1 — • — 1 — • — 

Pdec 

Pabs 

'"" defeated 


3.0 




2.5 



^2.0 
Q_ 

1.5 




1.0 




0.5 


1 1 AGeV (Au+Au) ■ 


3.5 


H 1 — 1 — • — 1 — • — 1 — • — 1 — • — 1 — • — 1 — • — 

Pdec 

- . Pabs 

: '-. defeated 


3.0 




2.5 
1- 
1.5 


"'"-C:: 

\ 


1.0 




0.5 


30 AGeV (Pb+Pb) ■ 


0.0 




3.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

M (GeV) 



FIG. 6: (Color online) Average baryon density experienced 
by a p meson as a function of the mass of the p meson. The 
results of the present calculation for the creation, absorption 
and decay point of the respective p meson are shown. From 
top to bottom we display central Au+Au/Pb+Pb collisions 
at 2, 11 and 30 AGeV. 



ever, this behaviour reverses already at i?i a b = HA GeV. 
At the highest energy (E^ = 30^4 GeV) displayed al- 
ready a factor of 2.5 more p mesons are produced in col- 
lisions than in decays of other resonances. For the loss 
term we observe a dominance of the absolute value of 
decayed p mesons at i?i a b = 30,4 GeV. 

Let us investigate the baryon density distribution at 
the space-time position of the p meson decay and ab- 
sorption in more detail. Fig. [4] shows the probability 
distribution of the baryon density at the instant of the 
p decay (dashed line), p absorption (dotted line) and for 
the sum of both (full line) . 

At 2A GeV (top figure) one clearly observes that 
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absorption of the p meson in the most dense medium 
reached at SIS energies (~ 3pb/pbo) is already a strong 
effect. If only decaying p mesons are taken into account 
the effective density probed reduces to (~ %Pb/ Pbq) At 
higher beam energies (11A GeV, 30A GeV) this splitting 
in the density between decaying p mesons and absorbed 
p mesons becomes even more pronounced. Here, absorp- 
tion processes are strongly dominating the high baryon 
density stage, while decays populate the low baryon den- 
sity region. 

To relate the present discussion directly to the sensi- 
tivity of di-lcptons to the most dense stages of the reac- 
tion, di-lcpton rates arc calculated as a function of the 
baryon density at which the di-leptons are emitted. The 
di-lepton calculations are based on standard cross sec- 
tions as e.g. given in [28|. Fig. Ogives the distribution 
of the baryon density at which the e + e _ -pairs from the 
p vector meson are emitted (from top to bottom cen- 
tral Au+Au/Pb+Pb collisions at 2, 11 and 30 AGcV are 
shown). For each energy, the di-lcpton production as a 
function of baryon density is separated into a part where 
the initial p meson decays (dashed lines) and a part where 
the p meson was absorbed by the medium (dotted lines), 
the sum of both contributions provides the total emission 
rate of di-leptons and is shown as full lines. 

Two different scenarios related to the initially dis- 
cussed theoretical approaches for the di-lepton extraction 
can be discussed: 

1. In the first scenario, the full branching ratio into 
di-lcptons is attached to both the p decay and the 
p absorption vertices (shown as thin lines). This 
provides the most optimistic reach towards high 
baryon densities, as it assumes that no collisional 
broadening takes place even in the most dense stage 
of the reaction (i.e. its a full vacuum baseline cal- 
culation). On the one hand, this setting is similar 
to the one employed in models that calculate di- 
lcptons from folding the y/s distribution of nucleon- 
nuclcon collisions as it assumes no interaction of the 
produced p meson with the medium, i.e. absorp- 
tion. On the other hand, if one omits the di-leptons 
emitted in the absorption process, these calcula- 
tions are similar to previous calculations that as- 
sume di-lcpton production only from the late stage 
decays of (baryon) resonances. 

2. This scenario is set in contrast to the shining ap- 
proach (indicated by thick lines). In the shin- 
ing method, a continuous emission of di-lcptons 
is assumed over the whole lifetime of the p meson 
(generally for all relevant hadrons). The di-lepton 
emission rate is then integrated over time, taking 
the collisional broadening for each individual vec- 
tor meson in its surrounding into account. Due to 
the strong collisional broadening in the medium, a 
drastic reduction in the analysis reach of di-lcptons 
towards high densities (shown in Fig. O from the 
comparison of the thick and thin full lines) results. 



At SIS and FAIR energies, the effective baryon den- 
sity probed by p mesons decaying into di-leptons 
is reduced to 1-2 pb/pbo m contrast to the ex- 
pected values around 2-3 pb/pbo from scenario 
1. At higher energies (11,4 GeV, 30,4 GeV), the 
reach of di-leptons into the most dense stage is also 
strongly reduced. In addition, late stage decays of 
baryon resonances and p mesons provide a strong 
trigger towards low baryon densities, resulting in 
strong peak of the di-lepton emission rate below 
1-1.5 pb/pbo- This low density peak might possi- 
bly blur the view on the most interesting di-leptons 
from the most dense stages of the reaction. 

Finally, we explore the average baryon density experi- 
enced by p mesons with different masses. This is impor- 
tant to understand whether any prominent features are 
present in the p meson mass region between 400 — 600 
MeV, that is important for the intermediate mass di- 
lcpton enhancement. Fig. [5] depicts the average baryon 
density experienced by a p meson as a function of the 
mass of the p meson. Fig. [5] gives the results of the 
present calculation for the creation, absorption and de- 
cay point of the respective p meson. One observes that 
the baryon density of the system is constant as a function 
of the mass, indicating that most p mesons, independent 
of their mass decay at a certain (low) baryon density, 
as argued before. Note that the baryon density where p 
mesons are absorbed is higher, which is in line with the 
previous discussion. 

In conclusion, we have shown that the measured di- 
lcptons provide only a restricted view into the most dense 
stages of the reaction despite the fact that electromag- 
netic probes leave the reaction zone without any further 
interaction. Thus, possible studies of meson and baryon 
properties at highest baryon densities might be blurred. 
For the p meson we have shown that the baryon den- 
sity probed in the di-lepton decay channel does depend 
on the method of di-lepton extraction employed. We ar- 
gued that the absorption of resonances in the high baryon 
density region of a heavy ion collision masks information 
from the early hot and dense stage. To demonstrate this, 
we have split the contributions of the loss term of the 
p meson yield into "loss(absorbed)" and "loss(dccayed)" 
and have shown that at early times, i.e. at the highest 
baryon densities the absorption results in substantial re- 
duction of the p meson life times. The explored energy 
regime of low and intermediate energies will be investi- 
gated with the planned critRHIC program and is in the 
focus of the future FAIR project. 
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Appendix: The four-current method 

The local baryon density at a space point i is the zeroth 
component of the baryon four-current — (pb,j)- The 
local rest frame (RF) baryon density at this space point 
is defined in the frame where the three-current vanishes, 
Jrf = (pb.rfJrf), with j RF = 0. This definition is 
known as the Eckart frame. Other definitions are possi- 
ble, e.g. in the Landau frame, the energy-momentum ten- 
sor is at rest while a baryon three-current might still be 
present. We believe however, that the Eckart frame defi- 
nition captures the relevant physics at the energy regime 
under investigation. 

In the context of the UrQMD model quantities are 
(per default) calculated in the computational frame (CF) 
which is (for symmetric systems) the center-of-mass 
frame of the whole heavy ion collision. In the com- 
putational frame one is only able to evaluate j^ F = 

{pb,cf,Jcf) where pb,cf = N/V is the baryon den- 
sity (N denoting baryon number in the volume, V be- 
ing the small local volume around the position i) and 
]cf = Pb,cfP- In the limit of an infinitely small vol- 
ume, the density pb,cf is a sum of Gaussians at position 
i: 

3 f (x-xo) 2 + (y-yo) 2 + (z-z ) 2 ~!z \ 

Pcf(t1) =Ef=i(7=^) 7,eV J 

V; ;/', ( 



i.e., a three-dimensional in ^-direction contracted and 
normalised Gaussian with j z = 1/ y/1 — /3| being the 
Lorentz factor for the particle under consideration. The 
normalisation is different for individual particles due to 
the different 7 factors. The nominal width of the Gaus- 
sian is case a = 1.5 fm. The particle that defines position 
i has to be included in the sum because one is interested 
in the baryon density in the local rest frame of the cell 
and not in the density around a particle in its rest frame. 

The velocity of the cell is computed with the same 
Gaussians as used for the density calculation as weight- 
ing functions. Therefore, the velocity of the cell in the 
computational frame is: 




The last step is to perform a general Lorentz boost of 
the four-vector into the local rest frame of the cell. 
I.e. a Lorentz transformation with the velocity of the cell 
@cf- The transformation matrix is the following: 

/ 7 -Pxl -Py7 -Pzl \ 

-& 7 i + ( 7 -i)f! ( 7 -i)^ (7-1)^ 

-A/7 (7-1)^ 1 + (7"1)$ (7-1)^ 

\-m (7-1)^ (7-1)^ i+(7-i)fiy 

(3) 

with (3 2 = Pi + $ + (3 2 z and 7 = - (3 2 . The zero- 

component of the transformed j M four-vector is the local 
) rest frame baryon density we are interested in. 
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